Intermittency in two-dimensional Ekman-Navier-Stokes turbulence 
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We study the statistics of the vorticity field in two-dimensional Navier-Stokes turbulence with 
a linear Ekman friction. We show that the small-scale vorticity fluctuations are intermittent, as 
conjectured by Nam et al. [Phys. Rev. Lett. 84 (2000) 5134]. The small-scale statistics of vorticity 
fluctuations coincides with the one of a passive scalar with finite lifetime transported by the velocity 
field itself. 



In many physical situations, the incompressible flow 
of a shallow layer of fluid can be described by the two- 
dimensional Navier-Stokes equations supplemented by a 
linear damping term which accounts for friction. An im- 
portant example comes from geophysical rotating flows 
subject to Ekman friction Q. The dynamics can be 
written in terms of a single scalar fleld, the vorticity 
a; = V X w, which obeys the equation 



'di 



(1) 



where v is the fluid viscosity, and a is the Ekman friction 
coeflticient. The term /^^ is an external source of energy 

- e.g. stirring - that counteracts the dissipation by vis- 
cosity and friction and allows to obtain a statistically 
steady state. Here, we will study the statistical proper- 
ties of vorticity fluctuations 5r^ = uj{x + r,t) — uj{x,t) 
at scales r smaller than the correlation length L of the 
external forcing. We will show that the statistics of 5ri^ 
is intermittent, and that the vorticity field has the same 
scaling properties as a passive scalar with a finite life- 
time. 

As shown in Figure ^ the vorticity field resulting from 
the numerical integration of Eq. (|l|) is characterized by 
filamental structures whose thickness can be as small 
as the smallest active lengthscales. The wide range of 
scales involved in the vorticity dynamics manifests itself 
in the appearance of power-law scaling for the spectrum 
of vorticity fluctuations Z{k) = 2nk{\uj{k)\^) ~ k^^~^. 
As already shown by Nam et al the spectral slope 

— 1 — ^ depends on the intensity of the Ekman drag: 
for the frictionless Navier-Stokes case {a = 0) we have 
^ = 0; a non- vanishing friction regularizes the flow de- 
pleting the formation of small-size structures and re- 
sults in a steeper spectrum (see Fig. ||). In the range 
< 1^ < 2 the exponent ^ coincides with the scaling ex- 
ponent C2 of the second-order moment of vorticity fluc- 
tuations S'^(r) = {{SrUj)^) ~ r^^. Let us now focus on 
a given value of a. In Figure |^ we show the probabil- 
ity density functions of vorticity fluctuations SrUJ at var- 
ious r, rescaled by their rms value {{Sr'-^yy^'^ . As the 
separation decreases, we observe that the probability of 
observing very weak or very intense vorticity excursions 



increases at the expense of fluctuations of average inten- 
sity. This phenomenon goes under the name of intermit- 
tency. Its visual counterpart is the organization of the 
field into "quiescent" areas (the patches, where vorticity 
changes smoothly) and "active" regions (the filaments, 
across which the vorticity experiences relatively strong 
excursions) . 




FIG. 1. Snapshot of the vorticity field. The numerical in- 
tegration of Eq. (0) is performed by a fully dealiased pseu- 
dospectral code with a second-order Runge-Kutta scheme, on 
a doubly periodic square domain of size L = 27r at differ- 
ent resolutions: A^^ = 512^,1024^,2048^ grid points. (The 
present figure is taken from the 1024^ run). Standard vis- 
cosity or hyperviscosity (depending on the resolution) is used 
to remove the remnant enstrophy flux at small scales. The 
large-scale forcing f^j is Gaussian, (5-correlated in time, and 
limited to a shell of wavenumbers around kf = 2-k /L {kf = 2 
for N = 512,2048, and fc/ = 8 for iV = 1024). At variance 
with other choices for fui commonly used (e.g. large-scale 
shear) this kind of forcing ensures the statistical isotropy and 
homogeneity of the vorticity field. 



The dynamical origin of this phenomenon can be un- 
derstood as follows (see also Ref. Q). Let us first notice 
that, for any a strictly positive and as far as the statisti- 
cal properties in the scaling range are concerned, we can 
disregard the viscous term in Eq. (|l|). In other words, this 
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system shows no dissipative anomaly, due to the presence 
of friction In the inviscid Umit {v = 0), Eq. (Q) can 
be solved by the method of characteristics yielding the 
expression uj{x, t) — /;j(X(s), s) exp[—a{t — s)] ds, 
where X{s) denotes the trajectory of a particle trans- 
ported by the flow, X{s) — v{X{s),s), ending at X{t) = 
X. The uniqueness of the trajectory X{s) in the limit 
^ is ensured by the fact that the velocity field is 
Lipschitz-continuous, as it can be seen from the velocity 
spectrum E{k) = Z{k)/k'^ ~ k~^~^, always steeper than 
k~^ (see Fig. ||). We remark that for ^ > the second- 
order velocity structure function is dominated by the IR 
contribution of the spectrum and thus trivially displays 
smooth scaling independently of the value of This is 
not the case for odd order structure functions that, in the 
absence of enstrophy dissipative anomaly, display anoma- 
lous scaling at the leading order Q. We have checked 
that this is indeed the case in our simulations. 



large times, the random variable 7 reaches a distribution 
P(7,t) ~ ii/2exp[-G'(7)i]. The Cramer function G{j) 
is concave, positive, with a quadratic minimum in A (the 
maximum Lyapunov exponent) G(A) — 0, and its shape 
far from the minimum depends on the details of the veloc- 
ity statistics Finite-time Lyapunov exponent and 
exit-times are related by the condition L — r exp[7Ti (r)] . 
That allows to obtain for r ^ L the following estimate 
for moments of vorticity fluctuations 



(2) 



The scaling exponents are evaluated from Eq. (g) 
by a steepest descent argument as = 
min-y {p, [pa -|- G(7)]/7]} . Intermittency manifests it- 
self in the nonlinear dependence of the exponents Cp on 
the order p. 
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FIG. 2. The vorticity spectrum Z{k) ~ k~^~^ steepens 
by increasing the Ekman coefficient a. Here a — 0.15 (-I-), 
a = 0.23 (x), Q = 0.30 (0). In the inset, the exponent ^ as a 
function of a. 
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FIG. 3. Probability density functions of normalized vortic- 
ity increments Srw/ {{SrUify^^ . Here, r = 0.20 (-I-), r = 0.07 
(x), r = 0.02 (v)- For large separations the statistics is close 
to Gaussian, becoming increasingly intermittent for smaller 
r. 



Vorticity differences are then associated to couples 
of particles ui{x',t) - Lj{x,t) = /*^[/c^(X'(s), s) - 
fuj{X{s), s)] exp[— a(^ — s)] ds. Inside the time integral, 
the difference between the value of /^^ at X' and that at 
X is negligibly small as long as the two particles lie at 
a distance smaller than L, the correlation length of the 
forcing; conversely, when the pair is at a distance larger 
than L, it approximates a Gaussian random variable Q. 
We then have SrUJ ~ f2 jljj'^^'^^ exp[— a(t — s)] ds ~ 
r2exp[— aTi(r)], where TL{r) is the time that a couple 
of particles at distance r at time t takes to reach a sep- 
aration L (backward in time). Large vorticity fluctua- 
tions arise from couples of particles with relatively short 
exit-times TL{r) <^ {TL{r)), whereas small fluctuations 
are associated to large ones. Since the velocity field is 
smooth, two-dimensional and incompressible, particles 
separate exponentially fast and their statistics can be de- 
scribed in terms of finite-time Lyapunov exponent 7. For 



It has to be noticed that the active nature of w has been 
completely ignored in the above arguments: the crucial 
hypothesis in the derivation of Eq. (H) is that the statis- 
tics of trajectories be independent of the forcing fi^. This 
is quite a nontrivial assumption, since it is clear that forc- 
ing may affect large-scale vorticity and thus influence ve- 
locity statistics, but it can be justified by the following ar- 
gument. The random variable ft arises from forcing con- 
tributions along the trajectories at times s < t — TL{r) , 
whereas the exit-time Tl is clearly determined by the 
evolution of the strain at times t — TL{r) < s < t. Since 
the correlation time of the strain is a~^, for TL{r) ^ 1/a 
we might expect that and T^ir) be statistically inde- 
pendent. This condition can be translated in terms of 
the finite-time Lyapunov exponent as r ^ Lexp(— 7/0;) 
and thus at sufhciently small scales it is reasonable to 
consider a; as a passive field. We remark that, were the 
velocity field non-smooth, the exit-times would be inde- 
pendent of r in the limit r — > and the above argument 
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would not be relevant. Therefore, the smoothness of the 
velocity field plays a central role in the equivalence of 
vorticity and passive scalar statistics for this system. 




FIG. 4. Snapshot of the passive scalar field, simultaneous 
to the vorticity field shown in Fig. |^ 

To directly check whether small-scale vorticity can be 
considered as passively advected by velocity, we also 
solved the equation of transport of passive scalar with 
a finite lifetime i,|7|,i 
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where the velocity field results from the parallel in- 
tegration of Eq. (^. The parameters appearing in 
Eqs. (0) and (||) are the same, yet the forcings f^^ 
and fg are independent processes with the same statis- 
tics. According to the picture drawn above, we expect 
to observe the same small-scale statistics for SrUi and 
Sr9^9{x + r,t)-e{x,t). 
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FIG. 5. Power spectra of passive scalar (x) and vorticity 
{+). Here a = 0.15. In the inset we show the same spectra 
compensated by k^'^''^ . 

In Figure ^ we show the power spectra of vorticity. 



Z(k), and of passive scalar Egik). The two curves are 
parallel at large /c, in agreement with the expectation 
C,2 — (2- We notice that the two spectra do not collapse 
exactly onto each other. At large scales we observe a big 
bump in Z{k) around k — kf which has not any corre- 
spondent in Eg{k). This deviation is due to the presence 
of an inverse energy flux in the Navier-Stokes equation, a 
phenomenon that has no equivalent in the passive scalar 
case. Due to this effect, the scaling quality of Sp{r) is 
poorer than the Sp{r) one, and a direct comparison of 
scaling exponents in physical space is even more difficult. 
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FIG. 6. Probability density functions of vorticity difi'er- 
ences {+) and of passive scalar ones (x), normalized by their 
respective standard deviation, at different scales r within the 
scaling range. 

However, we observe in Fig. ^ that the probability 
density functions of vorticity and passive scalar incre- 
ments, once rescaled by their root-mean-square fluctu- 
ation, collapse remarkably well onto each other. That 
proves, along with the result — C2 obtained from Fig||, 
the equality of scaling exponents of vorticity and passive 
scalar at any order: = Cp- 

The actual values can be directly extracted from the 
statistics of the passive scalar, which is not spoiled by 
large-scale objects. In Fig. ^ we plot the first exponents 

as obtained by looking at the local slopes of the struc- 
ture functions Sp{r). The numerical values for C,p are 
validated by the almost perfect agreement with the La- 
grangian exit-time statistics. 

In conclusion, we have shown that in presence of linear 
friction the small scale vorticity fluctuations in two di- 
mensional direct cascade are intermittent. Intermittency 
is a consequence of the the competition between exponen- 
tial separation of Lagrangian trajectories and exponential 
decay of fluctuations due to friction. Small-scale vortic- 
ity fluctuations behave statistically as a passive scalar, 
as it has been confirmed by a direct comparison. The 
smoothness of the velocity field appears to be a crucial 
ingredient for the equality of active and passive scalar 
statistics. 
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FIG. 7. The scaling exponents of the passive scalar (+). 
We also show the exponents obtained from the exit times 
statistics (0) according to {exp[—apTL{r)]) ^ r^p with aver- 
age over about 2 x 10^ couple of Lagrangian particles. The 
errorbars axe estimated by the r.m.s. fluctuation of the local 
slope. In the inset we plot the Cramer function G{'y) com- 
puted from finite time Lyapunov exponents (symbols) and 
exit time statistics (line). 

This work was supported by EU under the contracts 

HPRN-CT-2000-00162 and FMRX-CT-98-0175. Numer- 
ical simulations have been performed at IDRIS (projects 



011226 and 011411) and at CINECA (INFM Parallel 
Computing Initiative). 



[1] R. Salmon, Geophysical Fluid Dynamics, Oxford Univer- 
sity Press, New York, USA (1998). Other well known ex- 
amples are the Raylcigh friction in stratified fiuids, the 
Hartmann friction in MHD (J. Sommeria, J. Fluid Mech. 
170, 139 (1986)) and the friction induced by surrounding 
air in soap films (M. Rivera and X.L. Wu, Phys. Rev. Lett. 
85, 976 (2000)). 

[2] K. Nam, E. Ott, T.M. Antonsen, P.N. Guzdar, Phys. Rev. 
Lett. 84, 5134 (2000). 

[3] M. Chertkov, Phys. of Fluids, 10, 3017 (1998). 

[4] D. Bernard, Europhys. Lett. 50, 333 (2000). 

[5] E. Ott, Chaos in Dynamical Systems, Cambridge Univer- 
sity Press, Cambridge, UK (1993). 

[6] T. Bohr, M.H. Jensen, G. Paladin and A. Vulpiani, Dy- 
namical Systems Approach to Turbulence, Cambridge Uni- 
versity Press, Cambridge, UK (1998). 

[7] K. Nam, T.M. Antonsen, P.N. Guzdar, E. Ott, Phys. Rev. 
Lett. 83, 3426 (1999). 

[8] Z. Neufeld, C. Lopez, E. Hernandez-Garcia, T. Tel, Phys. 
Rev. E.,61, 3857 (2000) 



4 



